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Particles in a magnetic field and plasma analogies: 
doubly periodic boundary conditions 


P.J. Forrester 

Department of Mathematics and Statistics, 
University of Melbourne, Victoria 3010, Australia 


The iV-particle free fermion state for quantum particles in the plane subject to a perpendicular 
magnetic field, and with doubly periodic boundary conditions, is written in a product form. 
The absolute value of this is used to formulate an exactly solvable one-component plasma 
model, and further motivates the formulation of an exactly solvable two-species Coulomb gas. 
The large N expansion of the free energy of both these models exhibits the same 0(1) term. 
On the basis of a relationship to the Gaussian free field, this term is predicted to be universal 
for conductive Coulomb systems in doubly periodic boundary conditions. 


1 Introduction 


The two-dimensional one-component plasma is a model of two-dimensional point charges in equilibrium. 
The point charges all have the same sign and magnitude, +q say, and the system is neutralized by a 
uniform smeared out negative background, of total charge —qN. The total potential energy then consists 
of a particle-particle contribution, a particle-background contribution, and a background-background 
contribution. Using the fact the solution of the two-dimensional Poisson equation in free boundary 
conditions is given by — log |ri — r 2 1, for a system of N point charges with coordinates fq,..., Fat confined 
to a neutralizing disk of unit charge density centred about the origin, the corresponding Boltzmann factor 
is readily calculated to be proportional to (see e.g. ini) 

e-rEf=UUlU2 jFfc - r,f, T := q^ksT. (1.1) 

^<j<k<N 


Introducing the complex coordinates Zj = Xj iyj for rj = and recalling the Vandermonde 

determinant formula 


det[2^- 


i<j<fc<Ar 


one sees that jni in the case P = 2 can be written 


where 


Mz) ■■= 


( 1 . 2 ) 


(1.3) 

(1.4) 


Up to normalization, this is the absolute value squared of a free Fermi system of N particles for which 
the single particle wave functions in state k are proportional to HH). 

The wave functions 03 are realized as degenerate eigenfunctions in the lowest energy level (Landau 
level) for a single quantum particle confined to the xy-pleaie and subject to a perpendicular magnetic 
field B = Bz, B > 0. To revise this point (see e.g. ^), one notes that the corresponding Hamiltonian is 

H ^ ^{-ihS/+ -A)^ (1.5) 

2m c 


1 




where the vector potential A must satisfy 


V X A = Bz. 


The mechanism for the degeneracies in dini is the fact that H commutes with the quantum analogue of 
the square of the classical orbit centre, where 


X ^ Ily , 


y = y + 

n 


with I := \Jh,cleB the magnetic length and 

d 


IIj, = —ih 


dx 


By = -ih— 
dy 


■ 


generalised momenta. A standard calculation shows that the energy levels of are = (n + ^)hwc, 

Wc := eBimc, for n = 0,1, 2,... while the eigenvalues of X^ + F^ are equal to (2m + 1)^, m = 0,1, 2,- 

Moreover, for the ground state energy Eq = ^hwc the eigenfunctions of X^ + are 


ipmir) 


(27r;22™;2'"m!)^/^ 


( 1 . 6 ) 


The states csi are mutually orthogonal, and have the interpretation that they have definite values for 
the distance from the origin to the centre of their cyclotron orbit (implied by the eigenvalue of X^ + F^) 
which increases with m. The most dense X-particle state ip, in which the particles are fermions but 
otherwise non-interacting, is therefore obtained by constructing a Slater determinant from the states 
ipQ(r), ■ ■ ■ ,ipN-i{r), and thus up to normalization its absolute value squared corresponds precisely to to 
that specified by and dia. 

The aim of this paper is to develop the plasma analogy for the same quantum problem as above, but 
now with doubly periodic boundary conditions. This suggests a solvable version of the doubly periodic 
two-dimensional one-component plasma. For this particular model (which involves an X-body potential), 
the leading finite size correction to the bulk free energy can be computed exactly. Having a solvable doubly 
periodic plasma system identihed by the quantum problem, we consider the corresponding version of the 
two-dimensional Coulomb gas (mobile positive and negative point charges) at the coupling F = 2. The 
leading finite size correction to the pressure in the expansion of the grand potential can be computed 
exactly. It is found to coincide with the correction term obtained for the one-component plasma. We 
end with a discussion of the universality of this term for two-dimensional Coulomb systems in doubly 
periodic boundary conditions and in their conductive phase. 


2 The one-component plasma in doubly periodic boundary con¬ 
ditions 

The two-dimensional Poisson equation reads 

^ + ^ = -2-^(^--')%-y'). (2.1) 

Making use of complex coordinates, in free boundary conditions it has the solution 

^z,z') = -\og\z-z'\ (2.2) 
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We seek the solution satisfying doubly periodic boundary conditions 

^{{x + L,y),{x\y')) = ^{(x,y), {x\y')) 

^{{x,y + W),{x',y')) = ^{{x,y),{x\y')). 

For this task, following [7|, consider the Jacobi theta function 

OO 

9 i{z;q) = -i (-l)"9("-i/2)'e2*(n-i/2)z 

n— — oo 

OO 

= sin ^ n (1 “ (l " (l - . 

n=l 


(2.3) 

(2.4) 


(2.5) 


The fact that 9i is an entire function which vanishes if and only if z = ttto + nm, m,n G Z, q = e®' 
(Im(T) > 0), and that 9i{z; q) ^ z9{(0-q) as z ^ 0, tells us that 


$(z, z') := — log 


L\9i{Tr{z - z')/L;q)\ 


q := 


( 2 . 6 ) 


7r6»((0;g) 

satisfies the Poisson’s equation HD for 0 < x,x' < L, 0 < y,y' < W, with the further specification that 
112.211 holds as jz — z'j —>■ 0. Moreover, since 

0i(z + 7r;g) =-6'i(z;g) and 6»i(z + 7re; g) =-g“^e“^“0i(z; g) (2.7) 

we see that (ESI obeys (ESI, while must be modified to read 


TT 


$((a;, y + W),{x ,y)) = --{2y + W) + $((a;, y), (x', y')). 


( 2 . 8 ) 


This latter point is of no surprise, as it is not possible to solve the Poisson equation in doubly periodic 
boundary conditions unless it is made charge neutral and thus modified to read 


9^$ , , 27r 

From the above working, the doubly periodic solution of EiSIi is seen to be 

Hz,z') = ^ + Hz,z'). 


(2.9) 


( 2 . 10 ) 


From the viewpoint of the magnetic analogy, it is of interest to construct a 2dOCP from the quasi- 
periodic potential ESI) rather than the fully periodic potential iimni . Consider then a charged system 
of N mobile particles, charge +g, confined to the rectangle 0<x<L, 0<y< W, interacting via the 
pair potential EH). Also present is a smeared out uniform background of total charge density —N/LW. 
The corresponding particle-background potential is given by 


N „w 

U 2 ■= where V(z)= / dx' / dy'^{z,z') 

Jo Jo 


( 2 . 11 ) 


while for the background- background interaction we have 


U3 :=- 


jIJL [ 

2 LW Jo 


pW 

dx / dyV{z). 
Jo 


For the integral in EUl), according to ESI) we must evaluate 







Use of the product form in (E2Il shows 


i{y') = 



12 


With Ui denoting the particle-particle interaction, the Boltzmann factor is thus 


^-/3(Ui + U2+I73) 


l^ 7rd{(0;g) y^^^ 


Xg-(rArV6)log(fl;(0;9)/2)g-.prEf=ife-M//2)= Y[ \ 0 i{ 7 r{zk 




( 2 . 12 ) 


3 Magnetic analogy 


3.1 The A^-particle wave function 


We seek to construct the analogue of the states (ESI, and then to compare OS with the Boltzmann 
factor For the first of these tasks, we follow m The double periodicity associated with 

has as its fundamental domain a rectangle. In discussing the Hamiltonian EHl in doubly periodic 
boundary conditions it is more natural to take the fundamental domain as a parallelogram with corners 
at (0, 0), (0, L), {L + Wi, IU 2 ), (lUi, IF 2 ). But independent of this detail, it is not possible to construct a 
vector potential which is doubly periodic and satisfies V x H = Bz. Instead two vector potentials 
and , related by the gauge transformation A^ = A^ + V/, are defined. The vector potential A^ is 
periodic under x ^ x -\- while A^ is periodic under x 1 -^ x + Wi, y >—!■ y + II 2 . From these potentials 
corresponding ground state solutions ■0^ and can be constructed which satisfy 

tp^ix,y) = ijj^{x + L,y) and tp^{x,y) = {x + Wi,y + W 2 ), (3.1) 


and furthermore are related by 

Pj^{x,y) = 


(3.2) 


This latter fact follows from A^ and A^ being related by a gauge transformation. Note that the absolute 
values of and '0^ are equal and have the periodicity of the parallelogram. 

For the vector potentials we take 


= -Byx, = I ( - y) x + (a. - y) , 

which have the periodicity required by EU and are related by a gauge transformation with the scalar 
function given explicitly by 


fix,y) 






2 2 

X 


xy 


m 2 
-y 


2 v2iur ' 21U2' 

Substituting the property 13.211 in the periodicity equation for 0^ gives 

ip^{x,v) = ip^ix + Wi,y + 


(3.3) 


(3.4) 


But 0^ is periodic in x of period L so for solutions with the properties ED the magnetic held must be 
such that 


IU 2 


2ttPN 

L 


(3.5) 


for some = 1, 2,.... Since the parallelogram has area LW 2 , this condition says that the total magnetic 
hux BLW 2 is an integer multiple of the hux quanta 4>o = hcje. 
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An easy to establish general property of the Hamiltonian is that the lowest energy state periodic 
in X can be written in the form where f{u) is a Laurent series in u, and so we can 

write 

OO 

^ — 27rinz/L 


^ {x,y) = e~ 


E 


o„e 


(3.6) 


Substituting (TOni in (Toil allows a simple recurrence for an to be obtained. There are N independent 
solutions which when substituted back in 13.611 gives 




( 7 r(Tm - Nz/L)-q^) , 


where q := with t = {—Wi + iW 2 )/L, m = 0 , 1 ,..., A^ — 1 and 

OO 

e^(u-q):= Y. 


(3.7) 


(3.8) 


From the N single particle states (TTTIl . the A^-particle free Fermi state in US is formed. To relate 
it to the Boltzmann factor 12.1211 requires the following generalizations of the Vandermonde identity [5]. 

Proposition 3.1. With {q^]q^)oc = OiliCl 

fj,{q) := NN/2^-iN-l)iN-2)/24^^2. ^2yJN-l)iN-2) ,2 ^ 


For N odd 


det 




63 (y{xj + a - l/N)-,q^^^ 

N 

= 613(71^(0;^+a);g)/Ar(9) ei(y{xk - Xj)\q'^ 


(3.9) 


while for N even 


det 


1=1 


l<j<k<N 




_(y{xj +a- l/N);q^/^'j 

N 

+o^)'-:<i)fN{<i) n -Xj)]qy 


j—l ^^j<k<N 

In E3), JJT^) . 6»1 is specified by lg.51) , O 3 by i3.8\) . while 


e^{u-q)= Y (-!)>”' 


(3.10) 


(3.11) 


The determinants in Proposition Kh ll can be transformed into the form required by the free Fermi state 
in IE3, say, corresponding to E3. This is done by multiplying both sides of 13.911 with a = 0 by 


det[e2"/fe/7Vl ^ ^lV/2^(iV-l)(37V/2+l) 

and multiplying both sides of Jinn with a = — 7 rr /2 by 

det[e2-dfc+l/2)/7V] ^^iV/2.JV+l.(iV-l)(3JV/2+l)^ 

k=o .iV-1 

As a consequence the A^-particle state '0°^, defined as a determinant of theta functions 

1 


{{xi,yi),... ,{xN,yN)) ■■= 


VTr\{Li^f/^ 


N 


1=1 


xdet 


g-27ri(fc-i)2,/L0^ ( 7 r(T(fc - 1) - Nzj/L); q^) 




(3.12) 
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can be written in the factorized form 




i(^-i)(3^/2+i)/^(g) 


i=i 


where s = 3 for iV odd and s = 1 for N even. 


N 

x6»s( - Tr^^/L;^) 6»i (-7r(zfe - %)/L; g). 


(3.13) 


l<j<k<N 


We see that with Wi = 0, W 2 = W, 1/P = 2-k is closely related to the Boltzmann factor 

with r = 2. In fact the two expressions are proportional except that in (imil there is a factor of 
g27rWEj=i yj^ while in this factor is replaced by |0s(— Zj/L;q)\‘^. Equivalently, since 


= -27rpEf=iy) 


N 


i=i 


= p-2^Pi:/Li{y:i-Wl2f 


N 2 

0i{7rJ2i^j-iL-^W)/2)/L;q)\ (3.14) 


1=1 


we see that differs from 12.1211 with r = 2 by a constant factor times the many body term 


N 

- (i - iW)/2)/L-q'^ 


1=1 


(3.15) 


3.2 A doubly periodic plasma with A-body potential 


Being a free Fermi state, the many body quantum system corresponding to (linUl is exactly solvable in 
the sense that its /-point density matrix can be expressed explicitly as an I x I determinant. The diagonal 
term of this density matrix gives the /-point ground state correlation function. This correlation function 
is identical to that for the corresponding classical state with Boltzmann factor proportional to It 

follows that if the plasma system with Boltmann factor at r = 2 is augmented by multiplication 

by 13.1511 . a solvable model is obtained. Moreover, what was a quasi doubly periodic system now becomes 
fully doubly periodic. 

Rather than consider the correlation functions, let us consider the free energy of this system. The 
state ijj^^ is normalized so that integrated over 0 < < L, 0 < ?/j < IT {j = 1,N) gives 

unity. It thus follows from (ITT^ that 

pL pL pW pW 

/ dxi • • • / dxN / dyi - - dy^ 

^0 Jo Jo Jo 

X n 1^1 W^fe-^,)A;<z)l'=^!(iA(2p)-'/")'^(/^^(g))-^ (3.16) 

^<j<k<N 

Hence for the partition function of the solvable plasma system we have 




Zn 


^ pL pL pW pW 

mi, """X 'io 

y.^-0(Ui+U2+U3) 

log(y;(0;,)/2) (i^(2p)-l/2)^ ( 
00 

7r^(2p)-A2qi/6 Y[{1- q^^f, q := 


ei(iTY.^z,-{L-tW)l2)/L-q) 


6 














and consequently the total free energy is given by 


j3F ■= - \ogZN 


AT 

^ log p/ 27 r 2 + 2 log -'?''')) 


(3.17) 


One notes that the corresponding free energy per unit volume in the thermodynamic limit is precisely 
that known for the plasma system at F = 2 in a disk pp, and also there is no surface free energy in 
keeping with the system being doubly periodic. 

4 The two-dimensional Coulomb gas 

At the coupling F = 2 the two-dimensional Coulomb gas of equal numbers of positive and negative 
point charges in the grand canonical ensemble has been solved exactly in various boundary conditions 
[S1|^[TT 1 [TC 1 [T7|. Flowever these solvability properties do not carry through to the case of doubly periodic 
conditions with pair potential (ITTm I7|. Instead, the experience gained from the study of the one- 
component plasma above suggests that we consider a model with the quasi doubly periodic potential 
(ESI, and extend it to a doubly periodic model by multiplying by a term analogous to 13.1511 . The latter 


is chosen to be 



where Wj and Zj denote the complex coordinates of the positive and negative charges respectively. The 
total Boltzmann factor atV = 2, Wm 2 say, is thus given by 



ni<j<fc<jv ^i(^fc - Wj-; g)6>i(zfc - Zj]q) 

nffe=i^iK- 


(4.1) 


We seek to evaluate the grand partition function 



For this we require a classical determinant identity due to Frobenius Uni (see also El). 
Proposition 4.1. With F specified in 


N 


94:(^{wj - Zj) - a;q)F{wi, .. .,wn]Zi, ...,ZN]q) 
i=i 




(4.2) 


-6fia]q)9i{wj - g) J 


Using (Oil with 0 = 0 shows we can write 


Wn 2 = ( 04 (O;g))"det 


Oat [K{wj — Zk)]j,k=i,...,N 

[K{zj — Wk)\i,k=l,...,N Oat 


(4.3) 


where 


Kiw - z) := ^^i(0;g) 9fiw-z;q) ^ 
L9fiQ\q) 9i(w - z]q) 
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The structure (Oil is familiar in studies of the two-dimensional Coulomb gas at T = 2 [^. It allows the 
Fredholm theory of integral operators to be applied m which tells us that 


^2(0 = (04(0; q)f det(l + CK) = (04(0; q)f n(l + C^a) 


(4.4) 


where K is the integral operator such that the corresponding eigenvalues A and eigenfunctions ipi{x,y) 
and ip 2 (x,y) are specified by the coupled equations 


pW 

/ dx 2 / dy 2 ^ 2 {x 2 ,y 2 )K{zi - Z2) = \tpi{xi,yi) 

Jo Jo 

pL pW 

/ dx2 / dy2^i{x2,y2)K(zi-Z2) = \tp2{xi,yi). 

Jo Jo 


(4.5) 


As an aside, we note that because (Ot satisfies the Poisson equation EH, and further for 0 < 
a;i,a:2 < A, 0 < yi,y 2 < bF, K{zi — Z 2 ) is analytic except at zi = Z 2 , we have 


d 9i{z - z') 
dz 01 (z — z') 


tt5{x — x')5{y — y'). 


Hence applying djdzi to the first equation in 114.511 . and djdzx to the second equation, reduces the coupled 
differential equations 


7r^0'i(O;g) 

L04(O;g) 

7r^0((O;g) 

L04(O;g) 


^ 2 {x,y) 

^i{x,y) 


X-^-ipiix,y) 

X-^'tp 2 ix,y). 


Up to a scaling of the eigenvalue A, these coupled equations are themselves equivalent to the two- 
dimensional free particle Dirac equation. Indeed, such a relationship between the two-dimensional 
Coulomb gas at F = 2 and the two-dimensional Dirac equation is well known |^; the special feature 
here is that the components of the wave functions ipi and 'ip 2 must be doubly periodic. 


We now seek the explicit form of the eigenvalues. For this purpose we compute the Fourier expansion 


of 


9{x,y;q) 


04(7r(x -k iy)/L;q) 
6 i{Tr{x -k iy)/L\q) 


(4.6) 


by determining the coefficients g-niy, q) such that 


00 

g{x,y-q)= ^ 

n— — oo 


In fact these are known from our earlier study 0, where with q = e it was shown 

,^, 04 (O;g) I -W<y<0 

gn[y,q) - 2*0/^(O;g) 1 - g-(2n+i) | q-i^n+i) ^ 0<y<W. 

From this the eigenvalues in (14.511 can be determined. 


(4.7) 


Proposition 4.2. Let y := 7 r( 2 n -I- 1)/L and v := 2^1 Then for each n G Z the eigenvalues are 
determined by the roots of the equation 


cosh + u2)i/2) -1 = 0. 


(4.8) 







Proof. 


We seek eigenfunctions of the form 


Substituting these forms and H4.6fl . and using the orthogonality of shows that H4.5|l reduces 

to 


Aa„(?/') = j dyb„{y)gn{y' - y;q) 

>>bn{y') = J dyan{y)gn{y - y'-,q)- 

Substituting the explicit form (ICT for gn and differentiating shows that 


which are to be solved subject to the conditions 

aniW) 


A^(a„(j/)e-(2n+i)y/L) ^ 

A^(&„(2/)e--(2n+i)y/L) ^ 

bn{W) 


a„(0) 


= 1 , 


^n(O) 


= 1 . 


(4.9) 


(4.10) 


Solving the equations in (Ol) in terms of a linear combination of exponential functions, and determining 
the unspecified constants according to (OUll . gives the condition lICT . □ 

Substituting the eigenvalues as specified by (lOll in iimi we see that the product over the roots can 
be carried out according to the following general formula [TTlITfJi. 

Proposition 4.3. For f{z) an analytic function of z with zeros at z = ^j, j £ Z, and a product expansion 
of the form 


we have 


Thus 




IF 
f{-c) 

/(o) ■ 


2 

^(C) = (^4(0; g)) n( 


n—1 

2 


(04(O;g)) n 


cosh(W(^^ + (27r^)^)^/^) — 
cosh(W^) — 1 

^2W{u^ + i2TTCffF_2Wfi 


1 _ f.-Wiu^ + {2n(:f yF . 


We know from a previous study El that 

OO 

47^ E + (” - 1/2) Vi")'/" -in- 1/2)/l) 


n—1 


1 - 


-f3PL+ — 

L—*oo dE 


(4.11) 


where P denotes the renormalized pressure (a renormalization, by way of a cutoff in the sum over n 
in ilCTTi is needed because of the short distance singularity of the Boltzmann factor for a positive and 
negative charge at T = 2). Thus for W, L ^ oo 




2W(Ai^ + (27rC)^)^/^-2U^/x LW/3P_l/6 

^ e g ; 


g := e 


-ti-W/L 


9 









Furthermore, since m 


04iz; g) = n (1 - _ ^2n) 

n—1 


we see that 



Hence, for W and L large. 


2 ^ /I _ g-W{^J.^ + (27^Cy) 

n(- 

n—1 


2n1/2 


1 - 


\ 4 ~ 

-) - 11(1 - 9 '”)'- 


n—1 


- logS2(C) ~ -LW(3P + 2 log (gi/12 _ ^2n)^ ^ (4^12) 

n—1 

Notice that the finite size correction term is precisely the same as exhibited by JTTTIl for the corresponding 
one-component plasma system. 


5 Universality 

A continuum viewpoint of general two-dimensional Coulomb systems in their conductive phase (such 
systems may include arbitrary short range non-Coulomb potentials) reveals a relationship with Gaussian 
field theory [TUI [T31 ITC| . The latter, in the case of doubly periodic boundary conditions exhibits the 
constant term present in (ITT7I) and KWi . and as such predicts that this term will be present in the 
expansion of the free energy of a general conductive two-dimensional Coulomb system in doubly periodic 
boundary conditions. 

To expand on these points, let us begin by recalling some of the arguments from CHI The underlying 
assumption is that the universal features of the grand partition function Sc of a conductive Coulomb 
system are correctly accounted for by the continuum functional expression 

Sc = J'Dpexp(^-^ J J dfdr' p{r)G{r,r')p{r')Y (5.1) 

Here p{r) is the continuum charge density and G{r,r') is the 2d Coulomb potential, defined as the 
solution of the Poisson equation 123 subject to doubly periodic boundary conditions. The measure Vp 
normalized so that for G{r,f') = 1, Sc = 1. 

The expression is the continuum analogue of a multidimensional Gaussian integral. As such, 
taking into consideration the normalization of T>p, it has the evaluation 

Sc = (det[G(r,r')])-'/". 

But from the Poisson equation G{f,f') = (—so this can be written as 

=c = (de,(-Tv»))‘'t (5.2) 

On the other hand, the expression (—occurs in Gaussian field theory. The partition 
function for such a theory is defined by the functional integral 

Zc = y exp ( - ^ y drj • (5.3) 

Here it is required that (^(r) has doubly periodic boundary conditions so as to be consistent with the 
Coulomb system, and the normalization of 'D(p is chosen such that if —is replaced by unity, then 
Zc = 1. This is evaluated by diagonalizing — V^, which shows 
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Comparing (lOl and itCT shows 


log Sc = - log Zq. 


(5.5) 


For the two-dimensional Gaussian free field in doubly periodic boundary conditions, Zc as specified 
by (15.411 has been evaluated by Cardy |3]. Ignoring the zero eigenvalue (this has its origin in the need 
to regularise the Poisson equation HU according to (ESI) for a doubly periodic solution to exist) it was 
shown 

OO 

logZG = 2log(^q^/^^l[{l-q^^)). (5.6) 

n—1 

This is the 0(1) term in (15.171) and (in^ . as consistent with (1^ . and so establishes the stated univer¬ 
sality property. 
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